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Abstract. - A one-dimensional long-range model of classical rotators with an extended degree of 
complexity, as compared to paradigmatic long-range systems, is introduced and studied. Working 
at constant density, in the thermodynamic limit one can prove the statistical equivalence with 
the Hamiltonian Mean Field model (HMF) and a-HMF: a second order phase transition is indeed 
observed at the critical energy threshold ec = 0.75. Conversely, when the thermodynamic limit 
is performed at infinite density (while keeping the length of the hosting interval L constant), 
the critical energy Ec is modulated as a function of L. At low energy, a self-organized collective 
crystal phase is reported to emerge, which converges to a perfect crystal in the limit e —>■ 0. To 
analyze the phenomenon, the equilibrium one particle density function is analytically computed 
by maximizing the entropy. The transition and the associated critical energy between the gaseous 
and the crystal phase is computed. Molecular dynamics show that the crystal phase is apparently 
split into two distinct regimes, depending on the the energy per particle e. For small e, particles 
are exactly located on the lattice sites; above an energy threshold e*, particles can travel from 
one site to another. However, e* does not signal a phase transition but reflects the finite time 
of observation: the perfect crystal observed for e > 0 corresponds to a long lasting dynamical 
transient, whose life time increases when the e > 0 approaches zero. 


The study of long-range interacting systems experienced 
renewed interest, in the last decades, due to the in- 
reasing computational resources made available on mod- 
rn computers. Improved numerical simulations allowed 
0 enhance the statistics over previous investigations, en- 
bling for new ideas to be tested and existing theory to 
le challenged mm- In a long-range system the two body 
lotential decays with the distance r as r““, the expo¬ 
nent OL being smaller than the dimension of the embedding 
space. Many different systems fall in such a wide category. 
Among others, it is worth mentioning self-gravitating sys¬ 
tems and their applications to astrophysics and cosmol¬ 
ogy OH], charged plasma dynamics [5115] . hydrodynamics 
of two-dimensional vortices [Till], spin-wave interactions 
[^ITO] with their implications for the physics of lasers m- 
As largely testified by recent advances, long-range sys¬ 
tems display a rich zoology of peculiar behaviors, ranging 
from ensemble inequivalence to out-of equilibrium dynam¬ 


ics, and constitute an intriguing arena for novel develop¬ 
ments, of both applied and fundamental interests. 

Several long-range models exist in the literature, that 
are in principle suitable for elucidating key aspects of both 
equilibrium and out-of-equilibrium dynamics. Simplified 
toy models and paradigmatic case studies have been in 
particular proposed, which allow to progress in the anal¬ 
ysis by imposing a substantial reduction in the inherent 
complexity of the inspected physical problems. Working 
along these lines, it is instructive to elaborate on the ingre¬ 
dients that drive the spontaneous emergence of coherent 
self-organized structures, as seen in real experiments. As 
an example relevant for the topic addressed in this Let¬ 
ter, crystal-like structures have been reported to occur in 
long-range physical systems, as e.g. in dusty charged plas¬ 
mas PlfH). Coulomb-interacting cold atoms [13] or Bose- 
Einstein condensates m- The phenomenology of these 
latter states is not clearly understood: experiments seem 
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to suggest the presence of an intermediate phase during 
the melting of the solid structure, which displays addi¬ 
tional ordering at different scales m- When long-range 
interactions are active, the interplay between potential 
and kinetic energy terms give rise to correlations that do 
not decay rapidly with the distance. This latter effect 
may eventually translate into the emergence of non trivial 
thermodynamic transitions and peculiar symmetries. 

Available long range models often lack the necessary 
degree of phase space complexity to effectively reproduce 
the characteristics of such patterns. For instance, in the 
celebrated Hamiltonian Mean Field (HMF) model the 
spatial extension of the phase space is bounded in size, 
each particle’s position being mapped onto the circle 
0 G [0,27r[. Furthermore, the notion of inter-particles’ 
distance is not explicitly accounted for in the framework 
of the HMF formulation. Other long-range solvable 
systems, albeit spatially extended, may lack of transla¬ 
tional invariance, like the model mm, or encode for a 
peculiar lattice metric, as it is the case for the so called 
a-HMF model mm- This latter has a predetermined 
skeleton that is embedded in the model itself, and which 
reflects in the observed structures as a natural imprint. 

To bridge this gap and eventually describe the appear¬ 
ance of crystal-like long-range structures, we will here in¬ 
troduce and characterize a new model of coupled rota¬ 
tors, closely inspired to the aforementioned a-HMF and 
HMF models. The model that we are going to introduce 
will spontaneously freeze into a crystal-like phase, at suf¬ 
ficiently low energies. When the energy gets increased, 
the self-organized crystal melts, resulting in a disordered 
gaseous medium. The Letter is organized as follows: we 
first briefly recall the thermodynamic properties of both 
the HMF and a-HMF models. We then move to present¬ 
ing the extended model, focusing on its equilibrium ther¬ 
modynamical features, and demonstrating the existence 
of a crystal phase, via combined numerical and analytical 
means. 

The HMF model has imposed itself in recent years as a 
paradigmatic model for the study of long range interacting 
systems m- Indeed, it displays a rich variety of features, 
most notably related to its out-of-equilibrium dynamics. 
Moreover, computing its statistical properties via direct 
simulations of the microscopic dynamics is both simple 
and fast. The HMF model is characterized by the follow¬ 
ing Hamiltonian: 

N 2 1 ^ 

i—1 

where qi and pi stand for the canonically conjugated vari¬ 
ables of particle i, and N denotes the total number of 
particles. The HMF model displays, at equilibrium, a sec¬ 
ond order phase transition at critical energy Sc = 0.5, 
with a characteristic order parameter defined as the global 
magnetization M = ^ ^■(cos(gi), sin(qi)). The associ¬ 


ated phase space is bounded: the positions of the parti¬ 
cles can be mapped on the unitary circle, because of the 
periodicity of the potential and the mean-field nature of 
the interaction. In order to explicitly take into account 
the decay of the interaction potential as due to long-range 
couplings, an extension of the HMF model referred to as 
the a-HMF model [TTlfTO] has been proposed in the liter¬ 
ature. In this generalized model, the rotators occupy the 
sites of a regular lattice and the Hamiltonian writes: 



I - cosjqi - Qj) 


I* - Jl 


( 2 ) 


where a is a free parameter, which enables one to adjust 
the strength of the interaction, and ||i — jjl represents 
the distance between the lattice sites i and j. The 
normalization factor reads N = (2/7V)“ -|- 2 
and ensures extensiveness. For hnite sizes, in order to 
keep the system statistically invariant by translation 
along the lattice, it is convenient to impose periodic 
boundary condition. To accomplish this, it is customary 
to confine the lattice on a circle of radius {N — I)/27r 
and then consider in m the minimum distance on the 
circle between any pair of selected sites. For a < 1, the 
a-HMF is thermodynamically equivalent to the HMF 
model [iBHniiin], albeit it displays a lattice ordering 
enforced by the model in its microscopic formulation. 


Starting from these premises, and to eventually ex¬ 
plore the possibility of generating a non trivial crystal-like 
phase, which is not the mere byproduct of the underlying 
discrete spatial support, we here propose a straightfor¬ 
ward generalization of the aforementioned a-HMF model, 
that we will refer to as to the /3-HMF. Particles are no 
longer forced to sit on specific lattice sites, but can freely 
explore the hosting support. In practice, we substitute in 
Hamiltonian m the term ||i — j\\ by \\qi — qj\\, a replace¬ 
ment which does not produce divergences in the force field, 
as long as a < 1. Introducing a continuum distance pa¬ 
rameter, allows in turn to extend the phase space that 
particles can explore while, at the same time, removing 
the artificial constraint of a fixed lattice support. The 
/3-HMF model is in some respects reminiscent of the self- 
gravitating ring model [28], but, as opposed to this latter, 
it does not require the introduction of an artificial regu¬ 
larization of the potential at small scales. In this Letter, 
we will always deal with periodic boundary conditions on 
a circle of length L. Then, we will consistently define the 
distance between two particles with position qi and qj as 
dij = minjlgi - qj\,L - \qi - qj\} := \\qi - qj\\ namely, the 
shortest arc-length on the circle. Based on the above, the 
Hamiltonian of the /3-HMF can be cast in the form 


N 




i=l 


El 
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N 
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1 - COs(g, - qj) 


(3) 


When 0 < a < 1 the system is formally long-range: 
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this latter parameter can be tuned to control the strength 
of the long-range interaction. A{N, L) is a normalization 
constant which guarantees extensiveness, following the 
Kac prescription [53], and reads A = The 

above relation can be rigorously derived in the large scale 
limit and holds in general, as confirmed by numerical 
simulation. As a consistency argument, we remark that 
the expected normalization factor are recovered both 
in the HMF (a = 0) and a—HMF (qi = 27rfc -I- for 
L —>■ oo, G [0,27r[ and k G [0,L[) limits. As a further, 
preliminary comment, we note that as Qi are assigned to 
label the particles’ positions, in the interaction potential 
one should write in principle cos(z(qi — qj)), where z 
is a suitable constant that makes the argument of the 
cosine adimensional. Without loss of generality, we set 
z = 1 in the following. Moreover, in this Letter, we will 
assume L = 2nl, I G N, and introduce the rescaled density 
po = N/l to simplify the notation. 

We start our analysis by focusing on the low temper¬ 
ature regime. In this limit we expect to obtain minimal 
energy states. From the two-body interaction potential, 
one can imagine that the system will self-organize so to 
fill distinct available minima. It is however not evident 
that the process should eventually materialize in a crystal 
like structure, condensation of all particles in just one 
minimum being a priori also possible. Indeed, the further 
apart the particles sit, the wider the potential minimum 
appears. Hence, for a fixed energy amount, particles are 
more prone to fluctuate if they populate distant sites, 
so endowing a global repulsive entropic force. Since we 
assume particles to populate a bounded circular domain, 
we expect that the entropic repulsion will oppose conden¬ 
sation, so favoring in turn the spontaneous emergence of 
the crystal state, at low temperature. 

To test the adequacy of this interpretative scenario, we 
carried out a campaign of 7V-body simulations. These are 
performed by using a fifth-order symplectic scheme m, 
with a time-step 6t = 0.05. The equation of motion can 
be derived from Hamiltonian We wish to emphasize 
that the explicit presence of the inter-particle distance in 
@ heavily increases the required computational cost, as 
compared to the simple HMF case study, or even the a- 
HMF model [55]). For this reason, our analysis is limited 
to a maximum of iV = 10^ particles. The time needed to 
converge to equilibrium gets clearly larger, as the popula¬ 
tion size is increased. As an additional remark, we notice 
that for N = 2 the system is integrable. In this limit¬ 
ing case, one can analytically access the phase space of 
the system, and use the information to test the numerical 
implementation. It is indeed extremely useful to bench¬ 
mark numerical and analytical outcomes to avoid technical 
pitfalls which might originate when the subtle condition 
I ki ~ 11 0, is dynamically met. Furthermore both total 

energy and momentum are constants of motion, due to the 



Fig. 1: A few tracer particles positions in the q-space as function of 
time. The analysis of the recorded trajectories shows that a priori 
three dynamical regimes exist for the /3-HMF model 0, depending 
on the energy per particle e. At £ = e* ~ 0.1, the perfect crystal 
like structure is lost, and the particles can wonder to neighboring 
sites. Above ec ~ 0.75 the self-organized crystal melts, resulting in 
a disordered gaseous medium. Here N = 500, L = N and a = 0.5. 


translational invariance, and were numerically monitored 
over time. All simulations presented in this Letter assume 
an initial Gaussian distribution of momenta, with total 
momentum P = 0. For any given value of the energy, 
the system is observed to converge to the same equilib¬ 
rium state, independently of the distribution positions q, 
assigned at t = 0: we tested qi = 27rfc, qi = 0, qi from a 
uniform distribution, or Gaussian centered either g = 0 or 
on each q = 27rk, which was the fastest to converge. 

To characterize the thermodynamic fate of the system 
one can in principle rely on the magnetization M = |M|. 
When particles are uniformly spread over the hosting do¬ 
main, we have M = 0. Gonversely, M = 1 is equally asso¬ 
ciated to bunching on a single cluster, the condensed phase 
to which we alluded above, or to a distribution peaked 
on sites which are 27r-periodic, namely the sought crystal 
state. To remove this degenerancy and discriminate be¬ 
tween a crystal ordering phase and a condensate solution 
with all particles close to each other in space, one needs to 
compute the distribution of relative particles’ positions. 

We first start by qualitatively looking at some parti¬ 
cle trajectories, q versus time, t, as obtained through di¬ 
rect 7V-body simulations. Those are displayed in Fig. [T] 
Three rather different equilibrium regimes are apparently 
detected. At very low energy densities, the system self- 
organizes into a crystal in which particles are bounded in 
the vicinity of what appear to be lattice sites. When in¬ 
creasing the energy per particle above a threshold amount 
that we termed e* ~ 0.1, the system enters a partially 
disordered phase: particles are nevertheless localized over 
finite windows of time, and the global magnetization is 
still positive; a “soft” crystal phase emerges. At higher 
energies, for e > Sc 0.75, the system loses any traces of 
regularity and the crystal melts into a gaseous medium. 
As we shall prove in the following, only two equilibrium 
regimes exist: a self-organized crystal with no particle-site 
correspondence and a gaseous disordered phase. In other 
words, s* does not signal a genuine transition but indi- 
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rectly stems from the finite time of observation assumed 
in the simulations, as in the spirit of the Arrhenius law. 
The perfect crystal exists only in the limit e —)• 0, when 
the energy barrier between two adjacent site diverges: the 
perfect crystal displayed for 0 < £ < e* corresponds to a 
long lasting dynamical transient, whose life time increases 
when the zero energy limit is eventually approached. 

To cast this observation on solid grounds, we calculated 
analytically the single particle distribution as follows a 
maximization of the thermodynamic entropy, subject to 
the constraints of the dynamics. As compared to the HMF 
and a-HMF models, we have one more degree of freedom 
in choosing the thermodynamic limit. The system now 
depends on both the rescaled length I and the average 
particle density po- When performing the limit N ^ oo 
we can consider po oo at fixed I, or, conversely, the limit 
Z —>• oo at constant po- Indeed, the dependence on both 
Po or L can be made explicit in the scaling constant A. 
In the dual limit, when the density is kept constant and 
Z > oo, it is not so straightforward to analytically recover 
the analogue Hamiltonian. The problem lies in that the 
potential is modulated by the cosine with a periodicity of 
27r, but it is not periodic. So we are not allowed to perform 
trivial simplifications without any knowledge on the form 
of the distribution function. Our strategy will be to solve 
the system for finite length, and observe the scaling of the 
associated thermodynamic quantities for / —oo. 

To tackle the infinite density limit {N oo with 
I constant), we consider the one-particle density func¬ 
tion f{q,p,t) and denote the particle density as p{q,t) = 
f f{Pi <Zi t)dp. f{q,p, t)dqdp represents the fraction of par¬ 
ticles in the interval [q, q + dq\[p,p + dp] at time t. The 
system conserves the energy density 


^ = e[/] = ^ (y ^ 

where F represents the phase space volume, and 


V{qA) 



1 — a 
2 -kL 



1 — cos{q — q') 

\\q-dT 


dq . (5) 


Also note that due to translational invariance, the total 
momentum P[f] = /p f(q,p, t)pdqdp is conserved together 
with normalization condition Q[f] = f{q,pA)dqdp = 
1. We assume then a strictly Boltzmannian entropy 
S[f] = -/p/( 9,p,t)log(/(q,p,t))dgdp. To obtain the 
equilibrium function f{q,p) which maximizes the en¬ 
tropy S[f], under the constrains of the dynamics im¬ 
posed, we have to solve the maximization problem 
in.a^f{S[f]\J^f{q,p,t)dqdp = l,e[f] = E,P[f] = P}. By 
introducing the Lagrange multipliers /3, 7 and p, and re¬ 
quiring stationarity of the entropy functional, one readily 
obtains: 

( 6 ) 

where and D = e~^~^ and the equilibrium condition 
leaves V and / time independent. In the above calcu¬ 
lation we made use of 7 = 0, as P[/] = 0, for our 


class of initial conditions. Expression ([5]) is the standard 
Maxwellian distribution, where /3 is the inverse temper¬ 
ature and p stands for the chemical potential. As for 
the HMF model, one obtains an implicit system, since 
V{q) depends in turn on f{q,p). To solve the problem 
we use an iterative numerical scheme introduced in [28) . 
which generalizes the standard Newton technique by en¬ 
suring that entropy steadily increases towards the maxi¬ 
mum value at each iteration step. From equation (jb]) we 
find that the equilibrium density function p{q) becomes: 

p{q) = f f{q,p)dp = . 
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Fig. 2: Analytically predicted equilibrium distribution densities p(q) 
for different dynamical regimes and a = 0.5. Here we plot a small 
portion of the q axis with 5 periods for a distribution function com¬ 
puted for I = 100 periods. The crystal lattice appears as expected. 

A comparison between the theoretically predicted 
p{q) and the result of numerical simulations can be 
easily drawn for s < when particles occupy the 
self-emerging lattice, and for e > Sc, when particles 
are homogeneously spread. In both cases, we observe 
a good matching between theory and simulations, even 
with a low number of particles and low densities (data 
not shown). In the intermediate phase e* < e < ec, 
averaging effects due to particles relocations makes the 
direct comparison more cumbersome, since one has 
to cope with the the fact that due to translational 
invariance we have as well local vibrations of the crystal. 
To overcome these limitations, in figure |2] we confront 
the analytical predictions for the distribution of the 
inter-particle distances C(d) = f p{q + d)p{q)dq to 
the simulated data. C{d) is by definition translational 
invariant and the agreement between theory and simula¬ 
tions is almost perfect, for all the explored energy regimes. 

As an additional test, we proceed to compute the curve 
for the equilibrium magnetization M{e), for different val¬ 
ues of the scaling exponent a and for I = 1000. A second 
order phase transition between ferromagnetic M > 0 and 
paramagnetic M = 0 phases is found to occur at a critical 
energy density Ec, see Figure 01 Remarkably, for a < 0.6 
(for the chosen value of 1) the data collapse on the HMF 
transition curve and Sc = 0.75. For larger values of a the 
curves shifts towards the left and, consequently, Sc gets 
reduced. 

As anticipated the numerically determined energy 
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Fig. 3: The theoretical distribution C{d) of inter-particle distances 
for infinite density (solid lines) is compared to direct simulations for 
po = 1 (symbols). Here I = 500 and a = 0.5. Different curves corre¬ 
spond to distinct energy densities e. The horizontal axis is rescaled 
in a mod (27r) representation and the curves are normalized so to 
have unitary area. C(d) is translationally invariant: an almost per¬ 
fect matching is observed. 



Fig. 4: Comparison between the analytic magnetization curves and 
the data recorded via A'^-body simulations, for a random flat initial 
condition over the space positions. Data refers to Z = 1000. An 
average over a time window At = 5000 is performed. As expected 
data is noisy close to the critical transition energy Sc- The black 
thick curve represents the HMF equilibrium curve. 


threshold e* which putatively separates the perfect and 
soft crystals, does not leave a trace in the magnetization 
curves. For both large and small values of the scaling ex¬ 
ponent a the agreement between theory and simulations is 
convincing. A dedicated campaign of investigations sug¬ 
gests that, for all values of a, the transition energy Sc con¬ 
verges to the asymptotic value 0.75, in the limit ( —1 oo. In 
this limit, the transition curves of the /3-HMF and HMF 
(and clearly a-HMF) superpose to each other. In the fol¬ 
lowing we will provide an independent argument to show 
that the /3-HMF reduces to the a-HMF, when I oo. 

Let us start by simplifying the expression of the dis¬ 
tance term which appears in Hamiltonian ([3]). Define 
Qj = 2'Khj + Sj, with hj G {0, ••• ,n} playing the role 
of a lattice index and 0 < < 27r. Given the equilibrium 

distribution p{S, h) computed above, in the limit L -G oo, 
we ignore the contribution to the potential coming from 
pair of particles i and j which are close and have hi = hj. 
Under these working assumptions, for hi ^ hj the distance 
term in equation ([3]) can be simplified as: 

1 


^ ^ ~ ^j) 


2'Khi 

2'Khj 

OL Cl 

r 

2Khi 

2Khj 

L 

L 


L 

L 


The second term in Eq. © expression vanishes, when 
performing the sum in the limit L ^ oo. One loses 
therefore the dependence on <5. In other words, the 
numerator and denominator of the potential in equa¬ 
tion o become independent and the potential scales as 
— cos((5i — (5j))/(||/ii —Since the equilib¬ 
rium distribution function is homogeneous in the index 
h and since each interval (of length 27r) in which we have 
sub-divided the system, contains the same number of par¬ 
ticles, we can reorder the sum in the potential by avoiding 
to individually count the particles which share the same 
index hi. We then perform the replacement hi —> i. Apart 


for a multiplying constant depending on poi we recover 
the form of the a-HMF potential ([2). Since the steps are 
independent of a, we speculated that the a-HMF equi¬ 
librium solution, equivalent to the original HMF model, 
is also found within the /3— HMF setting, for any value 
0 < a < 1, provided I is sufficiently large. 

We now return to Fig. [3] to discuss again the long- 
lasting dynamical regime of perfect crystal. Particles are 
confined on lattice sites, as enforced by the fact that p{q) 
displays a periodic set of minima with p ~ 0. For high 
enough energies the minimum pmin of p{q) lifts clearly 
above zero: even though the crystal seems to preserve 
a well defined structure, this latter reflects a collective 
behaviour. Fluctuations instigate an effective mixing 
between adjacent lattice sites. If a truly transition be¬ 
tween perfect and “soft” crystal existed, the gap between 
consecutive sites should close Pmin = 0 for some value of 
the energy density, so resulting in an infinite barrier for 
the particles to cross as in the spirit of the Arrhenius law. 
Due to the observed analogy between the ,5-HMF and a- 
HMF models, the distribution of densities {q modulo 27r) 
obtained for the /3-HMF model superposes exactly (data 
not shown) to those obtained for the HMF. The interest 
of this conclusion are twofolds. On the one side, this 
observation implies that the universality class of the HMF 
extends to embrace non trivial models, as the /3-HMF 
is. Then, more importantly for the aforementioned issue, 
Pmin ~ exp(—/3M) for q = IT, meaning that Pmin formally 
vanishes only at zero temperature. Hence, as anticipated, 
only at zero energy density (zero temperature), one finds 
a perfect crystal lattice, particles being de facto trapped 
on specific lattice sites. For a non zero, but still small, en¬ 
ergy density, we end up with a collective crystal structure 
with particles that can in principle move from one site 
to another, while preserving a macroscopic crystal-like 
self-organization. The lattice emerges as a collective 
property of the distribution, and the trapping is effective 
for finite, though potentially large, times. The observed 
energy threshold e* separating perfect and “soft” crystal 
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regimes reflects therefore the finite time of observation 
and the finite size of the system rather then representing a 
genuine transition between two distinct equilibrium states. 

Summing up, we have here introduced and studied a 
novel long-range model, to shed further light on those spa¬ 
tial self-organized structures that can emerge when long- 
range coupling are at play. A continuous inter-particles 
distance term adds a degree of complexity with respect 
to the HMF and a-HMF models, while reproducing sim¬ 
ilar thermodynamics behaviours. In the extended phase 
space of the /3-HMF it is possible to observe a more rich 
phenomenology of equilibrium features. Specifically, we 
unravelled two different states, a self-organized effective 
crystal phase and a gaseous phase at high energies, which 
could not be seen in the reduced phase space of previ¬ 
ous long-range models. The crystal phases manifests as 
a generalized self-organized state, in which the local lat¬ 
tice emerge from the collective behaviour of the particles. 
A solid to gas transition manisfests at a critical energy 
which is sensitive to the scaling size L, and converges to 
0.75 in the limit L oo. In the crystal phase, a dis¬ 
tinction between “perfect” and collective crystal can be 
observed for finite sizes and finite time. Moreover we 
were able to analytically obtain the equilibrium distribu¬ 
tion function and apply it to describe the observed system 
behavior. Assuming the validity of the latter solution, the 
system can be formally reduced to the HMF in the limit 
I oo, thus reproducing the same equilibrium transition 
between magnetized and homogeneous states. When op¬ 
erating at constant I, but infinite density, the /3-HMF re¬ 
tuns a different equilibrium transition curve as compared 
to the HMF model. Further studies should aim at evalu¬ 
ating the interest of the crystal phase predicted within the 
realm /3—HMF model, in light of its potential applications 
to real long-range systems of experimental relevance. For 
instance the fact there is no direct correspondence between 
emerging crystal sites and effective particle position could 
be reminiscent of supersolid properties. Another open 
question has to do with the possible presence of quasi¬ 
stationary out-of-equilibrium states which were 

non investigated in this paper. 

* * * 
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